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л.  Н.  СРЕТЕНСКИЙ 

РАСПРОСТРАНЕНИЕ  ВОЛН  В  УПРУГОМ  ПОЛУПРОСТРАНСТВЕ 

ПРИ  ДВИЖЕНИИ  ПРИЛИВНОЙ  ВОЛНЫ  ПО  ПОВЕРХНОСТИ 

БАССЕЙНА  КРУГОВОЙ  ФОРМЫ 

§  1.  Допустим,  что  к  поверхности  упругого  полупространства  при- 
кладываются нормальные  напряжения,  являющиеся  периодическими  функ- 
циями временя.  Мы  будем  предполагать,  что  эти  напряжения  отличны 
от  нуля  лишь  внутри  некоторой  области,  имеющей  форму  круга.  Резуль- 
таты, получающиеся  при  рассмотрении  этой  задачи  теории  упругости, 
дадут  возможность  проследить  за  распространением  упругих  волн,  возни- 
кающих при  движении  приливной  волны,  если  мы  будем  считать,  что 
приложенные  нормальные  напряжения  равны  величинам  гидростатиче- 
ского   давления,    развивающегося    при    прохождении    приливных   волн. 

Предположим,  что  поверхность,  подверженная  давлению,  есть  плос- 
кость 2  =  0.  Возьмем  цилиндрическую  систему  координат,  взяв  в  качестве 
ее  оси  07,  —  вертикаль,  идущую  вверх,  и  располагая  координаты  г  и  <р 
в  плоскости  2  =  0. 

Упругое  смещение   V  можно  разложить  на  две  части,  полагая 


7=7'  — §га(1Ф, 
причем  скаляр  Ф  и  вектор  V'  будут  интегралами   уравнений: 

^!^  =  1аф 


д^у  1    * »- 

_  =  ___сиг1сиг17'; 
а1у  V'  =  0. 

Эта   общая  система   уравнений   теории   упругих  волн  может  быть  пере- 
писана в  следующем  виде,  используя  цилиндрические  координаты: 

;.2^  =  ДФ,  (1) 

х2-^  =  Д«'---^_,  (2) 

х2_  =  Дг;'-_  +  __.  (3) 
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х^^  =  Дш',  (4) 

|1  +  1«^  +  «^  +  21  =  о.  (5) 

В  ЭТИХ  уравнениях  величины  м',  г;',  го'  суть  проекции  смещения  V'  соот- 
ветственно на  радиус  —  вектор  г,  на  перпендикуляр  к  нему  и  на  ось  02,. 

Определим  решение  этой  системы  уравнений,  обладающее  данной 
частотой  а  и  удовлетворяющее  следующим  пограничным  условиям  при 
г  =  0: 


[В 


,     г':>а. 


Функция  Р{г,'^)  считается  известной;  символ  т  служит  для  обозна- 
чения компоненты  напряжения.  Эти  пограничные  условия  могут  быть 
шписаны  через  величины  Ф,  и' ,  р' ,  ю'  так: 

Гди/    ,    ди'        ^  д^Ф\        ^  /Эр'    ,     \  дъи'         2    д^Ф  \       ^      .^. 

^"  =  2,(^— ,-^)-ХАФ  =  |^  ^     ^^^  (7) 

При  решении  задачи  мы  должны  принять  еще  во  внимание  условия 
излучения  —  все  упругие  волны  должны  уходить  в  бесконечность. 

§  2.  Рассмотрим  сначала  уравнения  (2)  и  (3).  Найдем  частное  реше- 
ние этих  уравнений  следующего  вида: 

где  к  есть  любое  положительное  число,  а  5  есть  целое  положительное 
число.  Неизвестные  функции  V  {г)  и  7  (/•)  будут  удовлетворять  системе 
уравнений: 


здесь 


а^  =:  А;^  -|-  х2а2 


Общее  решение  этой  системы  уравнений,   обладающее  конечными  значе- 
ниями при  г  =  О,  записывается  так: 

С/  =  СЛ+1(аг)-Ь1)Л_1(аг), 

IV  =  С/,^!  (аг)  —  1)Л_1  (а/-); 

здесь  С   11   О   суть  две  произвольные  постоянные.  Отсюда  мы  получаем 
такие  частные  решения  уравнений  (2)  и  (3): 


12  Л.  Н,  СР  ЕТЕ  некий 


Возьмем  теперь  уравнение  (4).  Это  уравнение  обладает  частным 
решением  следующего  вида: 

ю'  =  ^е^(а'+8ф)+'с2  ^^  (а/.) ,  (9) 

где  А  есть  произвольная  постоянная. 

Обратимся  теперь  к  уравнению  (5).  Подставляя  в  это  уравнение 
вместо  и',  V' ,  ю'  их  значения  (8)  и  (9),  мы  находим,  что  это  уравнение 
будет  удовлетворяться,  если  постоянные  С,  О  связать  уравнением: 

а(С  — /))  +  ^Л=0.  (10) 

Возьмем,  наконец,  уравнение  (1).  Нетрудно  установить,  что  это  урав- 
нение будет  иметь  частное  решение  следующего  вида: 

•)  :;-      ф  =  5е^И+8ф)+л'г/^(аг),  (И) 

где   В   есть  произвольная  постоянная  величина,  а  к'  определяется  фор- 
мулой: 

^'2  =  ^2  —  а2/г2. 

Составим  теперь  для  е  =  О  значения  напряжений  т^г,  'Сгф,  '^22,  отве- 
чающие частным  решениям  (8),  (9),  (И).  Пользуясь  формулами  (6)  и 
(7)  для  этих  напряжений,  получаем,  отбрасывая  общий  показательный 
множитель: 

Тгг  =  |А   (у  сиА-\-кВ  —  (х.к'В]/а-1  (аг)  —(-^(х.А  —  кС  —  сск'В)  7а+1{<хг)    , 
т,,  =  111  и^  о^А^кП  —  оск'в)  Л_1  {<хг)  -^(1о:А  —  кС—  ^к'Ву,^^{хг  )1 , 

-22  =  -^г  [4}хЛА'а2  +  (а2  -^  А2)  (Х/г2а2  _  2^к'^)]  Л  (аг).  (12) 

Если  постоянные  А,  В,  С,  В  связать  зависимостями 

^^А+кВ  —  (хк'В  =  0,    ^о^А  —  кС  —  о^к'В  =  О,  (13) 

то    напряжения   т^г    и    т^ф    будут    равны    нулю  на  всей    плоскости  2  =  0. 
Решая  совместно  уравнения  (10)  и  (13)  относительно  А,  В,  В,  нахо- 
дим 

А  =  -^С,     В=^-^-^С,     В^-С. 

к  акк  ' 

Отсюда  частные   решения    (8),    (9),    (И)  и  формула  (12)  примут  вид: 
и'  =  —2С/^  (аг)  е^(а<+8ф)+'^2,  ^'  = ——  С1,  (аг)  е^(«»'+«')+*% 

ю'  =—'^  С1^  (аг)  е*(«"+8Ф)+Л2^     Ф  =  —  5^1±^СЛ  (аг)  е*И+«')+*'^    (14) 


где 


8  ^  (2^2  -1-  х2а2)2  _  ^кк'{к^  Н-  х2а2).  (15) 
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§  3.  Перейдем  теперь  к  решению  поставленной  задачи.  Нам  требуется 
найти  такие  функции  м',  г;',  м;',  Ф,  чтобы  при  2  =  0  напряжения  1гг, 
х,ф  были  равны  нулю,  а  напряжение  Тгг  определялось  формулой 

'"^10  ,     г>а. 

Мы  примем,  что  функция  /' (г,  ср)  имеет  следующий  частный  вид: 

При  решении  этой  задачи  будем  основываться  на  известной  формуле 
теории  функций  Бесселя: 


^  а  ^а  5  и^  (а^)  Л  (аг)  /  (6)  «/5  =  /  (г). 


Для  применения  этой  формулы  будем  рассматривать  константу  С, 
входящую  в  формулы  (14),  как  некоторую  функцию  параметров  а  и  ^, 
определяемую  так: 

С  =  Одля^>а;      С  =  ?^ -^^Л  (а^)    для  е<а. 

[I-  о 

Заменим  в  формулах  (14)  величину  С  этим  ее  значением  и  проин- 
тегрируем результаты  по  параметрам  а  и  ^  соответственно  в  пределах 
О,  оо  и  О,  а.  Тогда  решение  поставленной  задачи  будет  даваться  сле- 
дующими формулами: 


оо  а 

и'  =  -|-^  ^е*^;(аг)«га^&/з(а|)/(?)с/е-е^("'+8'), 

о  о 

оо  а 

г;'  =  —  ^'  ^  ^  е^^1,  (аг)  (1(х  ^  и,  (аС-)  /  (5)  а^  •  еЧр^+^'^), 

0  9 

оо  а 

м;'  =  —  1  ^  ?^  е"^!,  (аг)  с/а  ^  Ы,  (а$)  /  (&)  ^5  •  е^(-*+'^) , 
о  о 

оо  а 

Ф  =  --[  "^"'^+^'^  е*'^Л  (а/-)  ^а  \  ^^8  (ае)  /  (I)  с/^  •  е'('^'+«<Р). 


^ 


(16) 


; 


Величина  6",  входящая  в  эти  формулы,  является  функцией  параметра 
а  и  определяется  формулой  (15),  которую  можно  переписать  следующим 
образом,  принимая  во  внимание  значения  величин  к  и  А': 


.У  (а)  =  (2а2  —  х2з2)2  _  4а2  |Ла^  _  х^ з2  у  „2  _  /^г^г 


(17) 


Надо  отметить,  что  по  соглашению  ]/а2  —  х2а2^0  для  |а|>ха  и 
усх.^  —  П^::''  ДЛЯ  |а|>>/га;  для  других  значений  а  величины  корней  опре- 
деляются по  непрерывности  и  во  всем  дальнейшем  берутся  значения 
этих  корней  с  положительной  действительной  частью. 
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§  4,  Общие  формулы  предыдущего  параграфа  принимают  особенно 
простой  вид,  доступный  для  всестороннего  исследования,  если  нормаль- 
ные напряжения  обусловлены  прохождением  приливной  волны  полусу- 
точного периода.  Для  полусуточного  прилива  5  =  2,  а  высота  подня- 
тия воды  С  в  круглом  бассейне  равна  величине  статического  прилива, 
т.  е. 


Отсюда  получаем  для  Хг-  такое  выражение: 


где  К  есть  некоторая  постоянная  величина. 

Таким    образом,    для    волн,    образованных    полусуточным    приливом, 
мы  имеем: 


«  =  ^(тТ- 


Интеграл   по  переменному  I,  входящий  в  формулы  (16),  может   быть 
вычислен;  мы  имеем: 


^  ^/г  (а^)  /  (?)  ^^  ==  |-  ^-^3  («а), 
о 

Отсюда  формулы  (16)  могут  быть  переписаны  так: 

и'  = К\  ^Цт-  е'^  ^2  (аг)  /д  (аа)  с?а. 

о 

со 

^''  ""  ~11Г  ^  5  Т"  ^'''■^2  ^'^^^  "^3  (аа)  (1<х., 

о 

оо 

ю'  = К\  ^^  е^^2  (а/*)  /д  (оса)  с?а, 

о 

оо 

Ф  =  -^К\^  ^-^  е'^' ^2  (аг)  /з  (аа)  а<х. 


.  .^г(о<-г2<р) 


(18) 


Преобразуем  эти  формулы  к  новому  виду;  во-первых,  заменим  в  фо] 

1 
муле  для  м'/2(аг)    через    -2"  [А  (а'') — ^в  (°''*)].  ^  затем,  во  всех  формулах, 

увеличивая  область  интегрирования  добавлением  участка  от — эо  до  О, 
заменим,  как  это  обычно  делается,  функции  Бесселя  от  аргумента  осл 
функциями  Ганкеля  первого  рода,  имея  в  виду,  что 

Н1  (ре-"^)  =  е(Р+1)'^  Н^^^  (р), 
Н^1\ре'^)  =  е-(Р+»)'^  Н^^^  (р). 
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Мы  получим  тогда  вместо  формул  (18)    такие  формулы 

аК    Г     лкк' 


и'  =  -  ^   \  ~  е^^  /з  (аа)  Н'^'^  («г)  ^а, 

— оо 
со 


V    = 


со 

т'  =  -^\^  -^  е^^  /з  (аа)  Я^^  (аг)  с/а, 


Ф 


оо 

2^^    3 


а2  +  Л2 


(1) 


е^'^/зМ^У'М^а. 


,  ^г(в<+2ф) 


(19) 


Отметим,  что  интегрирование  во  всех  этих  формулах  ведется  по  дей- 
ствительной оси  от  —  оо  до  — г,  затем  по  полуокружности  радиуса  з, 
с  центром  в  начале  координат  и  расположенной  в  верхней  полупло- 
скости и  затем  по  действительной  оси  от  точки  з  до  оо.  Этот  путь  ин- 
тегрирования должен,  кроме  того,  обходить  точки  Рели,  которые  по- 
лучаем, приравнивая  нулю  функцию  «5".  Соответствующие  значения  а 
дадут  нам  поверхностные  волны,  представляющие  собою  собственные 
колебания  полупространства  с  данной  частотой  о  и  с  периодом  тс  по 
отношению  к  углу  ф.  Эти  волны,  которые  можно  считать  взятыми  с  про- 
извольной амплитудой,  мы  устраняем  из  рассмотрения. 

§  5.  Для  изучения  колебаний  полупространства,  описываемых  фор- 
мулами (19),  введем  вместо  гиг  новые  переменные  Я  и  со,  полагая 

2  =  — Лз1п(о,     г  =  /?со8со; 

эти  переменные  являются,  таким  образом,  полярными  координатами 
точки  на  плоскости  проходящей  через  ось  02,  и  образовывают,  кроме 
того,  вместе  с  углом  ср  полярную  систему  координат  в  пространстве. 

Возьмем  сначала  первую  формулу  из  системы  (19)  и  запишем  ее 
так: 

оо 

и^  =  —^{'^К  /з  (оа)  е-'^й^!^" Щ^'^  (ай соз <о)  с/а. 

— оо 

Если  мы  будем  считать,  что  радиус  К  есть  величина  большая,  а 
угол  со  принимает  значения  в  пределах  от  О  до  ^тг  —  (0(,,  где  соц  про- 
извольно малое  число,  то,  в  силу  указанного  выше  выбора  пути  ин- 
теграции, мы  можем  под  знаком  интеграла  заменить  функцию  Ганкеля 
ее  асимптотическим  выражением: 


Я;^'^(аЛсо8оз)=  У 


(« 


к  С08 


Ь^) 


илН  соз  <о 
Таким  путем  мы  находим  для  и'  следующее  выражение: 

2       _  оо  

и'  =  ^  е  ^~  ""  ^2 ^  ^5^^:^/з  («а)  е«^<«.")  б/а, 

(1  Г      т:К  С08  <о    3  Л^  о  \      / 


(20) 


16[  ;.;'..,       Э'  ;:-■  Л.  Н.  СР.ВТ.БНСКИ'Й 


где 


Аналогично: 


N  (а,  со)  =  а.1  соз  о)  —  81п  (о  |/а2  —  х^а^.  (21) 


21а  К         ~1 
е 


1     . 

тгг 


\хЕ  СОЙ  о 


^/^О?^^   5^Л(«а)а-<--.<г.,       (22) 


и,'  =  ^  <,     о      1/      „^         *;  -^  л  М  «""<'■">  «г»,  (23) 

— оо 
оо 

\/-17^ ^  77^  «^3  М  е^^^"'-")  ^а.  (24) 

г      т:^^  соя  о     3  >/^  а  ^'         ^  ^      '  ^      '' 


_  1 

2\1 


Ф-^'^е'^""' 


где 


М  (а,  со)  =  аг  соз  о)  —  81п  У (^  —  /г^а^.  (25) 


Для  получения  асимптотических  формул,  изображающих  величины 
и' ,  у' ,  ю' ,  Ф  при  больших  Н,  надо  исследовать  функции  Л/ (а,  со)  и 
Л' (а,  со),  как  функции  комплексного  переменного  а.  К  этому  мы  сейчас 
и  обратимся. 

§  6.  Рассмотрим  сначала  функцию  71/ (а,  со).  Эта  функция  имеет  две 
точки  ветвления: 

а^  =  а/г,        «2  =  —  ок. 

При  рассмотрении  интеграла  (24)  мы  берем  ту  ветвь  многозначной 
функции  ]/а2  —  а^/г^,  которая  при  действительных  а,  больших  чем  ак, 
по  абсолютной  величине,  обладает  действительными  и  положительными 
значениями.  Благодаря  этому  путь  интегрирования  в  интеграле  (24) 
должен  идти  по  нижнему  берегу  разреза  — оо,  — а/г  и  по  верхнему 
берегу  разреза  а/г,  сю.  Между  точками  —  а/г  и  а/г  функция  М  (а,  со)  имеет 
чисто  мнимые  значения,  обращаясь  в  нуль  при  а  =  а/гз1псо.  Производ- 
ная 71/а(а,  со)  равна  нулю    при  он!  — — а/гсозсо,  причем 

М  (а',  со)  =  —  га/г;     М"^:  (а',  со)  =     .    /  „ — . 

На  некоторой  линии  (Г),  проходящей  через  точку  а'  и  встречающей 
ось  абсцисс  под  углом  45°,  мнимая  часть  функции  М  (а,  со)  будет  со- 
хранять постоянное  значение  — а/г,  а  действительная  часть  будет  ме- 
няться от  О  до  — оо  в  обе  стороны  от  точки  а'. 

Для  получения  асимптотической  оценки  интеграла  (24)  мы  можем 
интеграцию  выполнять  по  кривой  (Г),  и  асимптотическая  оценка  этого 
интеграла  будет: 


— /?У» 


Ф  = 


г  {2Н^  со';^  ь)  —  х2) /,  (астЛ  соз  о)  е    '^°^^  Г        2о/«  аш' ы 


о2  К'аЛ  соз  со  [(^2— 2/12  соз  2^)2 -}-2Л^  31П  2о>  соз  со  К  х2— Л2  соз  21и] 


г     ^2о/.  «Ш'со^^^ 


Выполняя  здесь  интеграцию,  мы  находим  после  преобразований  сле- 
дующее окончательное  асимптотическое  выражение  функции  Ф,  с  вос- 
становленным множителем  е^(<''+2Ф);  .      ■. 
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^^.  1дсо.  (— а-соз^со  — 1)/з(астйсо8<о) 

ф  =  ±й±-  .  1_2 I - X 

^■>^^^^  (  Ь.^  \2         /гЗ  /         Л» 

(1  —  2-2-  созрел  1  +  2  "з  з1а  2о>соз  ю  |/   1 ^  со8«со 


— 1{аНН—Ы — 2Ф) 


X д .  (26) 

Надо  отметить,  что  при  получении  этой  формулы  следует  считать,  что 
угол  (О  не  равен  нулю,  а  более  точно  —  надо  предполагать  угол  со 
большим,  чем  наперед  заданное  малое  число  е.  Это  ограничение  вы- 
текает из  того,  что  при  (0  =  0  точка  а'  совпадает  с  точкой  ветвле- 
ния —  а/г. 

Рассмотрим  затем  формулы  (20),  (22)  и  (23).  Для  нахождения  асим- 
птотических выражений  величин,  определяемых  этими  формулами,  рас- 
смотрим функцию  уУ  (а,  о>)  в  зависимости  от  комплексного  перемен- 
ного а.  Эта  функция  имеет  две  точки  ветвления 

сс^  =  ах,      аз  =  —  зх, 

и  мы  разрежем  плоскость  переменного  а  по  линиям  —  оо,  (х.2  и  а^,  оо, 
проводя  притом  интеграцию  по  нижнему  берегу  первого  разреза  и  по 
верхнему  берегу  второго  разреза.  Производная  Л^^  (а,  со)  равна  нулю  при 
а'  =  —  ах  со8<о,  причем  N  (а',  со)  =  —  гах    и    ТУ^»  (а',  со)  = г-д —  . 

Применяя  метод  перевала,  мы  сможем  на  основе  этих  данных  о 
функции  N  {х,  со)  найти  асимптотическое  выражение  величин  и',  V' ,  го': 


и  = 


81112  (^    у     СОЗ^  0> ^  ^з{ааУ.  СОЗ  со)  ^_{(дасД_о<_2ф) 


/г2 
2аК1 


V'  = 


(ЛХСГ  /  Д2  л  ' 

соз^  2со  —  21  з1п  2а>  соз  |/   соз^  о -^ 

С032  со  -  ^  /з  (асТХ  соз  со)  |(^д_„,_2<р, 


т2 


^Ах-ст"  ,  /       „  к^  Н^ 


соз 


соз^  со  —  ^ 


гс  = 


аК1 


соз2  со -^  /з  (астх  соз  о)  ^_{(охН— о/-2ф) 


соз' 


2со  —  2/  зш  2со  соя  со  |/   соз^  со -^ 


К 


Составим  по  этим  формулам  величины  поперечных  смещений  по 
трем  взаимно  перпендикулярным  направлениям,  образуемым  увеличи- 
вающимися   /?,   со  и  ср;    обозначим    эти    смещения    соответственно    через 

мя,  Ра,  ги^. 

Мы  получим: 


/к^ 
—2 соз^  со  .  7з  (астх  соз  со) 

„«  =  и;     ^„  =  -^;;;^ ■ /,.  ^      X 

соз2  2со  +  2  зш  2со  соз  со  |/    "^  —  соз^  со 
-— г(охй— с.<— 2Ф) 

х- — ^ ;  (27; 

2     труды  Морского  гидрофизического  ин-та,  т.  VI 


I^^  л.    н.    сретенский 


что  же  касается  ю^,  то  его  можно  положить  равным  нулю,  так  как 
оно  изменяется  обратно    пропорционально    квадрату  Н.   Заметим,  кроме 

—2 соз^со  определяется  как  число  действитель- 
ное и  положительное  для   соз  (0  =  0. 

Таким  образом,    при    больших    Д    смещения  частиц  упругого  полу- 
пространства, обязанные    распространению    поперечных    волн,    будут  на-  ] 
блюдаться    лишь    в    плоскости    меридиана    частицы    и    притом  лишь  в 
направлении,  перпендикулярном  к  радиусу-вектору  частицы. 

Из  формулы  (27)  следует  далее,    что  для  углов  со,  удовлетворяюпщх 
неравенству: 

/12 
СОЗ^  6)  <"  — 5- 

колебание  частиц  полупространства,  обладаюш;их  одинаковыми  Л  и  ср, 
находятся  всегда  в  одной  фазе  (или  в  противофазе);  для  углов  же  о, 
удовлетворяющих  противоположному  неравенству,  фазы  колебаний  ча- 
стиц упругой  среды  будут,  при  данных  Л  и  ср,  зависеть  от  угла  о> 
и,  следовательно,  волновые  поверхности  приобретают  довольно  слож- 
ный вид. 

/"^2 
—2 — соз^со  знаме- 
натель формулы  (27)  не  обращается  в  нуль  ни  для  одного  значения  со. 
В  этом  можно  убедиться,  вернувшись  к  первоначальному  выражению 
(17)  функции  45",  равенство  нулю  которой  дает  уравнение  для  опреде- 
ления скорости  волн  Рэли  V.  Приравнивая  нулю  указанный  знамена- 
тель, мы  найдем: 

2  1 


СОЗ''  со  = 
но  7<  ' 


1 


к 

Рассмотрим    теперь    волны,    определяемые    потенциалом  продольных 
смещений  (26).  Мы  имеем  с  принятой  степенью  приближения: 


;^  С082  о>  —  1     -Ь  2  — ^  31П  2о>  СОЗ  о>  1/    ^——2  соз^  со 


X 


X 


— (а/1Й— о(— 2Ф) 


/^ 


'       ''«>^0;  ^1|^  =  0.  (28) 


к  соз  ыд(^  '  к  дс 

Отсюда  мы  видим,  что  при  распространении  потенциальной  волны 
будут  наблюдаться,  на  большом  расстоянии  от  центра  бассейна,  сме- 
щения лишь  по  радиусу  В..  При  данных  7?  и  сэ  колебание  всех  частиц 
находится  в  одной  фазе  (или  в  противофазе);  амплитуды  же  колебаний 
зависят  от  угла  со. 

Из  формул  (27)  и  (28)  можно  вывести  далее  общее  заключение,  что 
частицы  упругого  полупространства  описывают,  вдалеке  от  бассейна, 
эллипсы,  оси  которых  располагаются  по  радиусу-вектору  частицы  и  по 
перпендикуляру  к  нему,  проводимому  в  плоскости  меридиана,  содер- 
жатцего  рассматриваемую  частицу. 
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Полученные  нами  асимптотические  формулы    пригодны  для  вычисле- 

1 
НИН  соответствующих  величин  при  углах  со,  не  достигающих  у  тг;  кроме 

того,  значения  угла  ы,  близкие  к  О,  должны  быть  также  исключены 
для  возможности  применения  метода  перевала.  Таким  образом,  обо- 
значая через  (Ор  и  со^  два  произвольно  малых  угла,  мы  сможем  опре- 
делить область  приложимости  наших  асимптотических  формул  следу- 
ющими неравенствами: 

(О^  <^  со  <^  -^  7Г  — ■  С0(). 

Следовательно,  для  определения  смещения  точек  поверхности  полу- 
пространства наши  формулы  не  пригодны.  Соответствующие  формулы 
можно  получить,  рассматривая  формулы  (18)  для  2  —  О  и  г  большого. 

§  7.  Предыдущий  анализ  оставляет  в  стороне  наиболее  интересный 
случай,  отвечающий  точкам,  принадлежащим  к  свободной  поверхности 
полупространства.  В  этом  параграфе  мы  рассмотрим  этот  случай,  осно- 
вываясь на  некотором  видоизменении  метода  асимптотической  оценки 
интегралов,  содержащих  большой  параметр. 

Рассмотрим  интеграл 

оо 

[^/{'.)К,ма'с,  (29) 

6 

где  Кд  есть  символ  функции  Бесселя  порядка  5  и  второго  рода  от  мни- 
мого аргумента.  Накладывая  на  функцию  /  (т)  те  же  ограничения,  как 
и  при  получении  асимптотического  разложения    интеграла 


\  /  ('с)  е-^''  (1^ 


о 


при  больших  г,    можно    показать,  что    интеграл    (29)  может  быть  пред- 
ставлен следующим  асимптотическим    разложением: 


°°  оо 


где 


оо  т 
^^=^^^^з(^)с/^ 
О 

I  ат  суть  коэффициенты  разложения  функции  /  {х)  в  ряд  вида: 

оо         !!} 


-Т— 1 


Аналогичное    асимптотическое    разложение    может    быть    установлено  и 
т,ля  интеграла 


^/(т)Я1^\^хг)сгт; 


2* 
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МЫ  имеем; 


где 


—  00  со        (ч 


8, 
г 
тп=1  ,  А: 


— оо  т 

о 

§  8.  Возьмем  формулы    (19)    и    положим    в    них    2  =  0.    Рассмотрим 
сначала  выражение   функции  Ф;  мы  получим: 

оо 

ф  3=  _  ^  е^(а<+2Ф)  С  "'  +  ^'  /з  {<^г)  /з  (аа)  <^а.  (30) 

О 

Найдем    асимптотическое    разложение    этого    интеграла    для    больших 
значений  г.  Обозначая  интеграл  этой  формулы  через  ^{^),  положим: 

^{^)  =  ^^'Ч^)  +  ^^'Ч^),  (31) 

где 

со 

О 

со 

о 
причем 

^(а)=-^^^/з(аа). 

Возьмем  сначала  интеграл    ^^^^(г);  имея  ввиду  поведение  функции  Н^ 
для  чисто  мнимых  значений  аргумента  и  связь  этой  функции  с  функци- 
ей К^  мы  можем  интеграл  ^^^\^)  преобразовать  к  такому  виду: 


00 


функция  /^  (гЗ)  переменного  р  удовлетворяет  условиям  предыдущего 
параграфа,  благодаря  че^му  мы  можем  найти  асимптотическое  разложение 
^'-^^  (г)  для  больших  г.  Мы  имеем: 

отсюда  асимптотическое  разложение  интеграла  ^  '  (/•)  будет  таким: 

^^^\^)=^,.:^  +  ...  (32) 

Обратимся  к  рассмотрению   интеграла  (^^^^г).    Преобразуем    этот  ин 
теграл  к  новому  виду  изменением  пути  интегрирования. 
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Первоначальный  путь  интеграции  начинался  в  точке  а  =  О,  шел  до 
ТОЧКИ  а  = /га  и  переходил  затем  на  верхний  берег  разреза  (/га,  счэ),  опре- 
деляемого функцией  У"а2  —  1г^а^\  вдоль  этого  берега  разреза  рассматри- 
ваемый корень  имеет  положительные  значения;  начиная  от  точки  а  =  хс 
и  до  точки  а  =  оо  корень  у^а^  —  х^а^  имеет  также  положительные  зна- 
чения. Присутствие  под  знаком  интеграла  функции  /Г2^^(аг)  позволяет 
составить  новый  путь  интеграции  из  мнимой  полуоси  от  точки  О  до 
точки  — сх)/,  из  участка  (оо,  /га)  действительной  оси,  проходимого  по 
нижнему  берегу  каждого  из  разрезов  многозначных  функций  "[/а^ —  /г^о^ 
и  ]/а2  —  х^а^  и  из  участка  (/га,  сю)  той  же  оси  и  проходимого  по  верх- 
нему берегу  тех  же  разрезов. 

В  согласии  с  этим  будем  иметь: 

— оо  Ла  оо 

(?<^>  (г)  =1\   Р  (^8)  Н^^  (г,3г)  сг?  -Ь  5  +  ^  .  (33) 

о  со         На 

Изучим  каждый  из  этих  трех  интегралов. 

К  первому  интегралу  может  быть  приложена  формула  асимптотиче- 
ского разложения,  установленная   в  предыдущем  параграфе.  Мы  имеем: 

— оо 

г   ^  Р  Щ  Н?'  (г^г)  ^  =  -  ^-^-^  .  .  .  (34) 

о 

Возьмем  теперь  второй  интеграл  правой  части  формулы  (33) .  Мы  можем 
записать  этот  интеграл  так: 

^°  ""^  (2а2  —  х2о2)  7з  (аа)  НЧ\аг)  их 


г  1    г  (2а^  — х^о'')./з(аа)^/^'"(аг)йа 

^  ~  2    5  (2а2  —  х2а2)2  —  4а2  \^(х^  —  х%2  ]А^^23Гл%2 

>  оо  ' 

^  "г  (2а2-х%2)/з(аа)Я^2)(аг)сга 

2    3    (2а2  —  х2а2)2  +  41а2  Уу.^о^  —  а.^  Ка^  —  кЧ'^  '  ^      ' 


Третий  интеграл    формулы    (33)    может    быть  преобразован  к  такому 
виду: 

?         1  "г  (2а2  —  х%2) /з  (аа)  Я^^Чо"-)  <^а 

О  ~~  2    ^  (2а2  _  х2о2)2  —  4гаа  >^х2а2— а2  V а."^  —  к'^п'^ 
На  Ло  ^  ' 


4-11; 1^ 

2   ^    (2а2  — X 


(2а2  —  х2с12)  /з  (аа)  Д^^)  („^^  ^^ 


хо 


(2а2  —  х2ст2)2  —  4а2  Уа.^  —  у.Ъ^  У<х^  —  к^а^  ' 

ли 

Отсюда  имеем 

^"        'г  /г?    (2а2  —  х2(т2)  /з  (аа)  Я^^)  (^г)  «2  Кх2о2  —  «2  V  а2  —  /г2а2 

^    "^  3   "^  ^^  О  (2а2  —  х2ст2)4  +  16а*  {а.^  —  к^а'^)  (х2а2  —  а2)  ^'^• 

оо  Ла  Ло 

Для  больших  значений   г    можно    Н^^  (аг)  заменить  ее  асимптотическим 
выражением: 


22  Л.     Н.     с  Р  ЕТВ  ИСКИ  й 


Отсюда  для  больших  г  имеем: 


Г-4-?  =  4е  ~  1^  '''1/АС''     (2а--  -  Х^а2)  ^2  у-^2^2  __  ^2  ]/-^2  _  Д2^2  ^^  (д^)  ^„  I 

^  ^  )  г       7ГГ  3      (2а2  _  х%2)4  +  16а*  (а^  —  к^а'^)  {у.^а''  _  «2)     '  ^      ГТ  '     ? 

оо        Ла  Ло 

Найдем  асимптотическое  разложение  этого  интеграла.  С  этой  целью 
преобразуем  путь  интеграции  (Ла,  ха),  идущий  по  верхнему  берегу  раз- 
реза {ко,  сх))  корня  |/а2  —  к^а'^,  в  путь  интеграции,  идущий  по  второму 
листу  многозначной  функции  |/а2  —  Н^а'^  и  состоящий  из  двух  прямо- 
линейных путей:   (/га,  /га   — оог)  и  (ха  —  с\эг,  ха).   Мы  получим: 

1_[    ^      а_     С    1 


Мы  имеем: 


оо         1ю  ^     /|о  ха— оог^ 

/|0 — оог  .^ 

)  У  2гЗ  ■  .(7''(2/г2  — х2)3 


^-  тгг       .^~ 


Г  2гЗ  а^х* 


Отсюда  получим: 

5  +  $  =  ^  %^  {^ /з(а-)  ^-'  <""-  ^'>-  щ^.  /з  (аЬ)  .-.-}     (36) 

оо  /кг 

Принимая  теперь  во  внимание  формулы  (30)  —  (36),  мы  находим,  на- 
конец, следующее  асимптотическое  разложение  потенциала  Ф: 


Ф 


с  помощью  этого  выражения  для  Ф  мы  получаем  значения  смещений  при 
распространении  продольной  волны: 


аг  [Л  о2г2        \(2/г2  —  х2)з ""  'у.''  I  '1 


1  аФ  _  8аА'  Гх2-  Л2  ^      7з  (аха)       .(^,,...з^_„,)         НЧ,{а}ю)     -г(лаг-2Ф-а<+^-) 


Г    Эф  ^  а^        [  X*  '     (2Л^  — х2)3 

Из  общей  формулы  (19),  определяющей  Ф  при  любом  2,  можно  найти 
смеп^ение ^г—    для  2=0;  мы  имеем  для  больших  п 

дФ        АаК    1     1x2— Л2   г   /        ч    —  г  (хог— 2=Р— о/— -^п) 
—  тг- = т-^\ 1 — Уо(ака)е  ^        + 


/з  (аЛо)        —г  (/юг— 2Ф-о<—  -^  тт) 


4  (2/г2  —  х2)  ^ 
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Пользуясь  общими  формулами  (19),  можно  найти  смещения  и' ,  V' ,  ю' 
в  точках  ПЛОСКОСТИ  2  =  0  для  больших  значений  г. 

Соответствующие  асимптотические  формулы  были  найдены,  но  нашей 
просьбе,  С.  С.  Войтом,  которому  автор  весьма  признателен  за  оказан- 
ную помощь  в  работе. 

Мы  имеем: 


м'  =  — 


ги'  = 


(Л     о2г2\  X*  ^  ^  2  (2/га  —  х2)2 '^  Р 


С  помощью  всех  этих    формул     мы    находим    следующие  полные  значе- 
ния упругих  смещений  при  больших  г: 


2аКУ  у.^  —  Н^   Г        /гх2  г  /    /    ч     " 

I 


•Д  ^         г(о<+2Ф) 


V  =  0, 
^-— ^|4(2Д2_х2)./зИ^)^  ^        + 

Эти  смещения  складываются  из  смещений,  определяемых  проходящими 
продольными  и  поперечными  волнами;  с  возрастанием  г  они  быстро  убы- 
вают. 

Представлено  14.  X  1954 
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Зге-Ьепзк!!,   Ьеоп1с1  М1ко1ае- 
т±сЪ. 

Назргоз'Ьгапепхе  уо1п  у 
ирги^ош  ро1иргоз-Ьгапз-Ьуе . . . 
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